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Abstract This study focuses on the transfer of bound water 
and liquid water in wood. The moisture changes and distri¬ 
bution of six wood species (three softwoods and three hard¬ 
woods) were investigated in the longitudinal direction 
exposed to long-term moisture sorption in static environ¬ 
mental conditions. Most species used for the experiment 
reached an estimated maximum moisture content, which 
indicated that there might be no significant hysteresis in the 
capillary pressure curve due to air entrapment. The experi¬ 
mental data for the different samples were found to vary 
considerably. Using initial values obtained by the Boltzmann 
transformation, the Levenberg-Marquardt method was 
used to determine the moisture diffusivity from measured 
moisture content changes with time and moisture profiles. 
The validity was ascertained by comparing the numerical 
results with the corresponding experimental measurements. 
There was a point of discontinuity and an abrupt change 
in the slope of the diffusivity function around the fiber sa¬ 
turation point, which might slow the numerical solution 
process. 

Key words Inverse method • Sorption • Moisture transfer • 
Parameter estimation • Finite volume method 


Introduction 

Wood is a typical hygroscopic porous material and consists 
of cell walls, cavities, and pits between the cell walls. Mois¬ 
ture in wood exists as three phases; namely, bound water, 
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water vapor, and free (or liquid) water. Wood is also sus¬ 
ceptible to growth of bacteria, fungi, or algae, especially 
when exposed to high humidity or liquid water. Liquid 
water may provide a medium for diffusion of enzymes or 
other metabolites by which biological deterioration can be 
facilitated. The deterioration of decayed wood may be 
accelerated by physical and chemical damage. The durabil¬ 
ity of wood in building largely depends on moisture condi¬ 
tions. For traditional post and beam wooden structures, 
moisture content is kept below the fiber saturation point 
(FSP) during general use. However, wooden columns may 
be directly exposed to rainfall and come into direct contact 
with liquid water. In addition, moisture may enter struc¬ 
tures by condensation, runoff from roof and facade, and by 
capillary action of groundwater. Therefore, simulation of 
moisture transport in wood under set environmental condi¬ 
tions should be valuable for predicting susceptibility to bio¬ 
logical decay. 

Liquid water uptake of porous materials is called imbibi¬ 
tion or absorption in hydrology and building physics. In the 
case of hygroscopic porous materials such as wood, however, 
the more general term might be sorption. Because the mois¬ 
ture content of exterior wood is in an air-dried condition, 
moisture changes in wood in contact with liquid water 
involves transport of bound water, water vapor, and liquid 
water during sorption. As a consequence, most models use 
Fick's law to model bound water and water vapor for the 
hygroscopic region that is below the FSP. Darcy’s law is 
used for free-water transport for the nonhygroscopic region 
that is over the FSP. The heat and moisture transfer coeffi¬ 
cients of wood in the hygroscopic range, under the FSP, 
were reviewed by Kang et al. u To model moisture sorption 
including in the over-hygroscopic range, liquid transport 
parameters such as capillary pressure and intrinsic and rela¬ 
tive permeability are required. The experimental data are 
very limited . 2 A considerable number of experimental and 
theoretical studies have been conducted on moisture trans¬ 
port in wood drying. However, most such models contain 
at least one adjustable parameter, which must be deter¬ 
mined by an experiment. The relevant heat and mass trans¬ 
fer coefficients are necessary for the simulation whose 
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accuracy depends mainly on the input coefficients. Informa¬ 
tion relating to moisture transport properties is limited and 
insufficient, especially for the simulation of the sorption 
process in wood. 

Moisture diffusivity can be estimated from moisture 
changes with time and transient moisture content profiles 
at a given time by inverse methods: (1) the Boltzmann 
transformation method (BTM), (2) the moisture gradient 
method (MGM), and (3) the optimization method (OM). 
By using the BTM, Koponen 3 investigated diffusion coeffi¬ 
cients of wood in drying over both hygroscopic and over- 
hygroscopic ranges under isothermal conditions. Diffusion 
coefficients increased with increasing moisture content up 
to the FSP and then decreased in the nonhygroscopic region. 
By using the MGM and the gradient in water potential as 
the driving force, Cloutier and Fortin 4 investigated the 
moisture transport coefficient of wood during drying under 
nonisothermal conditions in which it increased exponen¬ 
tially with moisture content. Hukka 5 calculated the mois¬ 
ture diffusion coefficient below the FSP from the measured 
development of the internal moisture profiles during drying 
by the conjugate gradient method of the OM. Due to hys¬ 
teresis between drying and sorption, these results are readily 
applicable to the wood sorption process. However, there 
has been little investigation concerning the measurement 
and modeling of moisture sorption in wood in which the 
moisture content increased from low to high. 6-9 

Descamps 10 reported that the optimization method with 
a spline function performs better than the BTM or MGM 
because it does not use the derivative of a moisture profile. 
However, the optimization methods have some serious limi¬ 
tations and different initial settings can lead to different 
optimization results. Sometimes the functions become over¬ 
parameterized and the same curve can be generated by 
different parameter combinations. It is also helpful to use 
expressions with parameters that can be directly obtained 
from the experiments. The goodness-of-fit depends on a 
good guess of initial parameters. 

The purpose of this study was to investigate the moisture 
change and distribution of six wood species (three soft¬ 
woods and three hardwoods) in the longitudinal direction 
when exposed to long-term moisture sorption. Based on a 
transport model of sorption, diffusivities of bound water 
and free water covering the whole moisture range needed 
in the simulation of wood is determined by using the inverse 
method. In this study, the initial settings obtained by the 
Boltzmann transformation method are used and compared 
for accuracy. Validity is ascertained by comparing the 
numerical results with the corresponding experimental 
measurements. 


Theory 

Moisture transfer model 

Multiphase heat and mass transfer models that considered 
moisture content, air, and energy conservation were devel¬ 
oped by several researchers. If the gas pressure in wood is 


assumed to be constant, the mass transfer model can be 
converted to the diffusion equation: 

^ = V{(D w + D^D^)Vm] (1) 


where m is the moisture content based on the dry weight of 
wood, D„ is the liquid water diffusivity, A is the water 
vapor diffusivity, and D B is the bound water diffusivity. 

The liquid water diffusivity (D w ) over the hygroscopic 
range can be expressed by using the generalized Darcy’s 
law: 


A 


Pw Aw-^w dPc 
p 0 V w dS dm 


Vm > m FSP 


( 2 ) 


where p 0 is dry wood density, p w is liquid water density, k IW 
is relative water permeability, K w is intrinsic water perme¬ 
ability, v w is dynamic viscosity of liquid water, P c is capillary 
pressure, and S is liquid water saturation. 

The combined diffusion coefficient (Ac) under the 
isothermal condition is expressed in different equations 
depending on the models. 2 


r. -73- TT~ _ (7 n V0(V0 -p)( l+1 -V^)Pc£>bl 

Ac - A- + Al - Pv 0 A +---/=— - — — — 77 ;- 

1 - V 0 + Lp c D BL /A, 

+ (1-0)PcAjl m<m FSP (3) 


where p is the ratio of an effective pit opening area to end 
wall area; L is the ratio of fiber length to an effective diam¬ 
eter, D y is the diffusivity of water vapor in air, D B is the 
diffusivity of bound water in the cell wall, <j> is the porosity 
or void fraction, and p c is the cell wall density. 

The porosity of wood may be calculated from wood and 
cell wall density using Eq. 4, assuming that the density of 
bound water is equal to the water density, the oven-dry cell 
wall density is 1540 kg/m 3 , and the pore volume remains 
constant with the moisture changes. 11 


p 0 1000 x Go 

1540 1540 


( 4 ) 


The bound water diffusion coefficient in the longitudinal 
direction c an be estimated from transverse one, 

D BT = 7 x 10" 6 exp [- (38500 - 29000m)/R T ] (5) 

Al = 2.5 x At (6) 

and the water vapor diffusion coefficient is 


IX = IX, 


0.018p vs dh 
p'"RT dm 


Av =2.2x10" 


f 1 . 013 -loA 

f T ) 

l p J 

V 273 J 


1.75 


(A 

(8) 


where R is the universal gas constant (8.314 J/mol K), T is 
temperature in Kelvin, p„ is saturated vapor pressure, h is 
relative humidity, and P is total pressure of air and water 
vapor. 

According to Siau, 11 Eq. 3 may be expressed in simplified 
form. 
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D bc = ——-- j- (9) 

1-0 d hl + (\-^)d v /l 

The maximum possible moisture content of wood can be 
predicted theoretically if the cell cavities remain constant 
in size during moisture sorption and m FSP is known. 


: '^FSP + 


( 10 ) 


where the fiber saturation point of wood wz FSP is usually 
assumed to be 0.3 and G 0 is the dry specific gravity of 
wood. 

However, it is difficult to measure the dry volume because 
the moisture content can change during measurement. 
Therefore, the dry specific gravity can be obtained using the 
moist one. 


1 -mG„ 


( 11 ) 


where G m is the specific gravity of wood based on dry weight 
and moist volume at moisture content m. 


Numerical analysis 


The equation governing isothermal moisture transfer can 
be obtained by substituting Eq. 3 into Eq. 1. 


dm 

Ik 


= V-[(D v/ + D bc )Vm] 
= V[Z)Vm] 


( 12 ) 


Finite volume method 


Consider the nonlinear diffusion equation as 
dm d 


dx . 


, t> 0, 0 <x<L 


dt dx . 
with the initial condition 
m(x, 0) = m 0 (x), 0 <x<L 
and two point boundary conditions 
dm 


D 


dx 


= h[m(L, t)-0.05] 


x=L 


m (0, t ) = m n 


(13) 


(14) 


(15) 


where h is the external convective mass transfer coefficient 

(0.0001/p o ). 

In order to find the solution of Eq. 13 by the finite volume 
method, one takes some meshes 

n = {(*;, f ; -)|0 = t 0 < q<... < t N , 0 = x 0 < Xi<... < x M = L} 

(16) 

and constructs an appropriate control volume correspond¬ 
ing to spatial nodes jx,j. 

By integrating Eq. 13 on each control volume (Fig. 1) 
and using the Gauss-divergence theorem, the following dis¬ 
crete form is obtained: 


-G- 1 0 

Xu Xj 

Control Volume 


-©— 

X i+ 1 


Fig. 1. Nodes and control volume 
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At 


X i D { ml+1 ) 

f=E,W 


3m ;+1 , 

dx J 


■n f 


(17) 


where superscript j is the time layer, subscript i represents 
nodes, and is the outward normal vector at each end point 
of the control volume. At the end points W, the value of 
D(m i+1 ) is given by an average value 


D{mi + \= {D{m\ +1 )+ D{m\tl))l2 


(18) 


and the derivative is approximated numerically as in Eq. 19 


dm’ +1 _ mj+l - mj +1 
dx x i+ \ - x ; 


(19) 


At the ( 7 +l)-th time, Eq. 19 provides the solution m <;) at the 
(y')-th time. Hence, as shown in Eq. 17, each equation can 
be used to obtain m[l\, m( +1 , and m’~l. However, the system 
of equations must be nonlinear, and so Newton’s iteration 
method has been used. 


Diffusion coefficients 

To adopt the generalized Darcy’s law to the simulation, 
there is a need for an intrinsic and relative permeability as 
well as a capillary-liquid water relation. For the combined 
diffusion coefficient, the law requires information of the cell 
wall geometry as well as the sorption isotherm, for which 
there are no sufficient data until now. Therefore, two types 
of diffusion coefficient were assumed as follows: 

First, one coefficient was proposed by the exponential 
function with the fourth-order polynomial, which is contin¬ 
uous in the whole moisture content range. 

D = eycp(a Q + a l m + a 2 m 2 + a 3 m i ^ (20) 

Second, two coefficients were proposed by the exponen¬ 
tial function with a third-degree polynomial over the hygro¬ 
scopic range and a first-order polynomial in the hygroscopic 
range. 

D w = exp(a 0 + a \ m + a 2 17,2 + a 3 m3 ) Vm > zzzfsp 

D bc = exp(-£> 0 -h 1 m) Vm<m FSP (21) 


Least-square method and Gauss-Newton method 

Suppose that one must determine the value of the coeffi¬ 
cient D(m) by using the known experimental results of m[ 
for some i and j. Because there are four parameters a„, a,, 
« 2 . and a 3 , many experimental values for m\ are given. In 
this case, it is common to use the least-squares method as 
follows: 
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To find the parameters A = [a 0 , a u a 2 , a 3 ] that minimize 
the objective functional, it was defined as the sum of the 
squares of the deviation between the measured and the 
predicted data on both average moisture content and a 
moisture profile. 


N 2 M l 

F(A) = m e {tj)\ + X [m(x k , t N )-m e {x k , t N )f 

/=1 k =1 

( 22 ) 

where in and m e are the average moisture contents of pre¬ 
dicted and experimental data with time t , respectively. 
m e (x k , t N ) is the moisture profile of the experimental data 
for final time t N and m(x k , t N ) is one of the numerical solu¬ 
tions for the nonlinear diffusion equation (Eq. 13) by finite 
volume method. 

The condition necessary to minimize the function at 
A* is 


VF(A*) = 0 

Then the gradient of F(A) is 
VF(A) = f> ; (A) Vr t (A) = J(A) T r( A) 

i =1 

where r(A) = (r u r 2 ,..., r K ) T 

r i (A) = in(t i )-ih e (t i ) i = l...N 

r N+k( A ) = m{x k ,t N )-m e {x k , t N ) k = l...M 
K = N + M 

and the Jacobian matrix of r(A) is 


f 

dr l 

dri 

I 

da 0 


da 2 

da 3 

dr 2 

dr 2 

dr 2 

dr 2 

da 0 

dfl! 

da 2 

da 2 

dr K 

dr K 

dr K 

dr K 

V da 0 

dcii 

da 2 

da 2 . 


(23) 


(24) 


(25) 


However, this system of equations is nonlinear and difficult 
to compute. Hence Newton’s method is used to solve this 
system. The original Newton’s iterative formula has the 
form 

A (* + 1) = A^~[j(A^) T J{A^)+ r(A«) 

(26) 

where Hessian S(A) has the second-order derivative such 
as 


S(A) = f j r i {x)V 1 r i (x) 

i =1 


(27) 


However, S(A) is expensive to compute and makes the 
system ill-conditioned. Hence, by neglecting this second- 
order term in Newton’s method, the simplified iteration is: 


A (k+l) = A (k) - 


j(A™) T J(A^)] 1 J(A^) T r(A«) (28) 


This equation is said to be the Gauss-Newton iteration. To 
make the iteration well defined, it is required that the rank 
of the Jacobian matrix is 4. Obviously, the success of the 
Gauss-Newton method will depend on the importance of 
the neglected second-order term S(A). If S(A*) = 0, the 
Gauss-Newton method would have a quadratic conver¬ 
gence rate. In general, the Levenberg-Marquardt method is 
used to improve the numerical stability and convergence. 
These parameters were adjusted by the fitting procedure 
until the best agreement between experiment and simula¬ 
tion was obtained. 


Materials and methods 

Three species of each softwood and hardwood were selected 
for the experiments (Table 1). Wood samples were cut to a 
prismatic shape with approximate dimensions of 25 mm 
(width) x 25 mm (depth) x 200 mm (height). Four vertical 
sides were coated by a varnish with heavy pigments. The 


Table 1 . Physical properties of wood species used for the experiment 


Wood species 

Duration 

(months) 

Dry density 
(kg/m 3 ) 

Porosity 

Maximum MC (kg/kg) 

zlMC a 



Estimated 

Measured 


Softwood 

Hemlock ( Tsuga heterophylla) 

1 

485 

0.69 

1.71 

1.80 

-0.09 


2 

448 

0.71 

1.88 

2.02 

-0.14 

Radiata pine ( Pinus radiata) 

1 

422 

0.73 

2.02 

1.94 

+0.06 


2 

461 

0.70 

1.82 

1.74 

+0.08 

Spruce (Piceci mariana) 

1 

510 

0.67 

1.61 

1.61 

0.00 


2 

500 

0.68 

1.65 

1.69 

-0.04 

Hardwood 

Painted maple (Acer mono) 

1 

626 

0.59 

1.25 

1.21 

+0.04 


2 

650 

0.58 

1.19 

1.03 

+0.16 

Japanese elm ( Ulmus davidiana) 

1 

800 

0.48 

0.90 

0.86 

+0.04 


2 

686 

0.55 

1.11 

1.03 

+0.08 

Horn beam (Carpinus cordata) 

1 

741 

0.52 

1.00 

0.90 

+0.10 


2 

724 

0.53 

1.03 

0.93 

+0.10 


MC, Moisture content 

“Differential moisture content between estimated and measured maximum moisture content 
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top and bottom faces were left uncoated. The bottom 
surface was submerged only a few millimeters in order to 
avoid buildup of hydrostatic pressure. The top surface was 
exposed to ambient conditions. The environmental labora¬ 
tory conditions were kept constant at 20°C and 40% relative 
humidity during the experiments (Fig. 2). 

For all samples, moisture content (MC) was recorded 
at various time intervals. The moisture of each sample in 
grams was obtained by subtracting the moisture of the 
oven-dried samples from the moisture of the wet samples. 
After the absorption test, the specimen was cut consecu¬ 
tively into 24 to 27 pieces for measuring moisture and spe¬ 
cific gravity distributions. The volume for the specific gravity 
was measured by the water immersion method after both 
uncoated surfaces were coated with paraffin. 


Results and discussion 

Maximum moisture content 

The estimation of maximum MC is important in predicting 
the sorption process of wood accurately, which should be 


Fig. 2. Schematic diagram of 
the moisture absorption test. 
MC, Moisture content 


MC 5.0% 



defined at the boundary conditions. Most maximum MC 
values estimated by using Eq. 4 were slightly higher than 
the measured values; hemlock was the exception. This may 
be due to lower FSP, air entrapment, or dead pores. The 
difference in hardwood was higher than in softwood. 
Flowever, Eq. 4 can be used to gauge an approximate value 
of maximum MC. 


Estimation of moisture diffusivity 

Moisture content changes with time, and moisture distribu¬ 
tions at the final stage are shown in Fig. 3. The moisture 
diffusivity values in Tables 2 and 3 were obtained by the 
Boltzmann transformation method over the whole moisture 
content. 12 Flowever, the input values of b Q and b l were 
approximated using Eq. 9. 

The shape of the curve representing the dependence of 
diffusion coefficients with moisture content was similar in 
all wood, and the diffusivity decreased with moisture content 
until the FSP was reached and then increased. Adapting 
them to the moisture transfer equation (Eq. 12), the simula¬ 
tion results gave poor predictions for moisture content 
changes with time and moisture distributions. This is because 
the Boltzmann transformation method is valid for a semi¬ 
infinite volume and a linear correlation between moisture 
flux and the square root of time. 

Mathematically, inverse problems belong to the class of 
ill-posed problems. The solution does not satisfy the general 
requirements of existence, uniqueness, and stability under 
small changes to the input data. 13 The inverse solutions are 
known to be very sensitive to small changes in input data 
resulting from measurement and modeling errors. There¬ 
fore, to improve the numerical convergence, the initial 
input values of parameters obtained by the Boltzmann 
transformation method were used as a first approximation 
for the Levenberg-Marquardt method. Equation 9 for the 
moisture diffusivity below the FSP was used initial input 
values for the discontinuous model. 


Table 2. Diffusivities with the continuous function by the inverse method, Eq. 20 


Species 

Duration 

(months) 

«0 

“i 

a 2 

a 3 

F 







LM“ 

B b 

Softwood 

Hemlock 

1 

-22.35 

3.571 

-7.875 

3.849 

0.161 

0.777 


2 

-22.37 

1.880 

-4.050 

1.969 

0.102 

0.160 

Radiata pine 

1 

-21.05 

-5.596 

8.226 

-2.714 

0.524 

0.608 


2 

-21.52 

-6.264 

8.237 

-2.276 

0.575 

3.001 

Spruce 

1 

-21.80 

-4.354 

8.995 

-3.552 

0.266 

0.377 


2 

-21.49 

-6.456 

8.889 

-2.592 

0.408 

0.438 

Hardwood 

Painted maple 

1 

-16.35 

-17.63 

24.09 

-8.392 

0.140 

15.53 


2 

-18.14 

-14.84 

28.25 

-13.37 

0.125 

7.092 

Japanese elm 

1 

-21.01 

-4.527 

7.630 

0.381 

0.146 

0.146 


2 

-21.62 

-0.857 

-6.190 

8.955 

0.361 

0.469 

Horn beam 

1 

-21.52 

-3.563 

-1.538 

9.700 

0.046 

0.246 


2 

-21.45 

-6.995 

20.57 

-12.78 

0.016 

0.020 


a Levenberg-Marquardt method 
b Boltzmann transformation method 
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Fig. 3a-f. Comparison of 
measured and predicted average 
moisture content change and 
moisture distribution using two 
steps of moisture diffusivity 
by the Levenberg-Marquardt 
inverse method, a Hemlock, b 
radiata pine, c spruce, d painted 
maple, e Japanese elm, f horn 
beam. Measured: filled circles , 1 
month; open circles , 2 months. 
Predicted: solid line, 1 month; 
dotted line, 2 months 


Distance from wetted surface (mm) Distance from wetted surface (mm) 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Nondimensional time Nondimensional time 


(a)Hemlock (b)Radiata pine 


Distance from wetted surface (mm) 

0 50 100 150 200 



0.0 0.2 0.4 0.6 0.8 1.0 

Nondimensional time 


(c)Spruce 


Distance from wetted surface (mm) 

0 50 100 150 200 



0.0 0.2 0.4 0.6 0.8 1.0 

Nondimensional time 


(d)Painted maple 


Distance from wetted surface (mm) Distance from wetted surface (mm) 

0 50 100 150 200 o 50 100 150 200 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Nondimensional time Nondimensional time 


(e) Japanese elm (f) Horn beam 


Table 3. Diffusivities with two step functions by the inverse method, Eq. 21 


Species 


Duration a 0 a 1 a 2 a 3 b 0 b t F 

(months) - 

LM B a 


Softwood 


Hemlock 

1 

-24.83 

9.567 

-11.99 

4.706 

17.95 

14.87 

0.042 

0.320 


2 

-22.58 

2.552 

-5.179 

2.407 

19.95 

10.35 

0.056 

10.22 

Radiata pine 

1 

-20.32 

-14.49 

18.70 

-5.887 

18.64 

12.66 

0.272 

2.055 


2 

-21.59 

-2.453 

1.298 

0.610 

19.65 

3.574 

0.333 

4.873 

Spruce 

1 

-23.76 

4.817 

-1.515 

-0.014 

17.68 

9.609 

0.118 

0.192 

Hardwood 

2 

-20.27 

-11.93 

14.85 

-4.384 

15.96 

21.40 

0.059 

0.664 

Painted maple 

1 

-19.55 

-10.29 

18.67 

-7.152 

17.72 

9.629 

0.114 

1.475 


2 

-25.12 

16.79 

-19.86 

10.34 

18.87 

3.010 

0.078 

8.431 

Japanese elm 

1 

-24.09 

17.41 

-32.35 

22.32 

13.74 

34.61 

0.064 

0.115 


2 

-23.18 

14.72 

-33.27 

22.41 

15.69 

19.33 

0.078 

0.901 

Horn beam 

1 

-29.98 

44.12 

-82.70 

52.86 

19.26 

11.15 

0.028 

0.153 


2 

-22.93 

2.491 

3.353 

-3.284 

19.37 

10.86 

0.037 

0.330 


“Equation 9 below the fiber saturation point and Boltzmann transformation method in nonhygroscopic range 
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Fig. 4a-f. Apparent moisture 
diffusivity over whole moisture 
content predicted by the 
Levenberg-Marquardt inverse 
method, a Hemlock, b radiata 
pine, c spruce, d painted 
maple, e Japanese elm, f horn 
beam. FSP, Fiber saturation 
point 








Moisture content (kg/kg) 
(f) Horn beam 


Comparing Tables 2 and 3, the discontinuous diffusivity 
model (Eq. 21) gave better prediction of the development 
of the internal moisture distributions in wood during sorp¬ 
tion than the continuous model (Eq. 20). Using the discon¬ 
tinuous diffusivity of Table 3, the simulation results are 
represented by solid lines in Fig. 3. 

This study shows that there is a distinctive hygroscopic 
zone in wood, although both experimental and model-based 
determination of phase-separated transport properties is 
still one of the most disputed issues in moisture transport 
analysis. 14 For liquid water sorption simulation, moisture 
diffusivity should be separated into hygroscopic and nonhy- 
groscopic zones and two step functions of diffusivity should 
be used rather than an effective diffusivity to improve the 
simulation results. 

As shown in Fig. 4, the moisture diffusivity decreased 
with increasing moisture content until the FSP was reached 


and it increased in the nonhygroscopic range, which agrees 
with other studies. 71115 For most species, the order of 
magnitude of moisture diffusivity is lies within the range of 
1CT 8 to lO" 10 m 2 /s. 


Conclusions 

Using the Levenberg-Marquardt inverse technique, two 
diffusion models to predict moisture change with time and 
moisture distribution have been investigated for liquid 
water sorption. The moisture-dependent diffusivity was 
considered as a continuous function over the whole mois¬ 
ture range, while this coefficient was discontinuous at the 
fiber saturation point. The initial input values of parameters 
by the Boltzmann transformation method were used to 
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improve the numerical convergence as a first approximation 
for the Levenberg-Marquardt method. Comparing the cal¬ 
culated moisture changes and profiles with experimental 
values, it is apparent that the discontinuous model using the 
parameters obtained by the inverse method is able to predict 
the development of the internal moisture distributions in 
wood during liquid water sorption. The results of this study 
should be useful in studying the process of sorption for 
exterior wood such as columns and dipping methods used 
for preservation treatments. 
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